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ABSTRACT 

^ ■ An analysis of the kinematics and ionisation state of the emission hne gas of a sample 

of 14 3CR radio galaxies with redshifts 2; ~ 1 is carried out. The data used for these 
studies, deep long-slit spectroscopic exposures from the William Herschel Telescope, 
0> ' are presented in an accompanying paper. It is found that radio sources with small 

linear sizes (;< 150 kpc) have lower ionisation states, higher emission line fluxes and 
broader line widths than larger radio sources. An analysis of the low redshift sample 
of Baum et al. demonstrates that radio galaxies at low redshift show similar evolution 
in their velocity structures and emission line ratios from small to large radio sources. 

The emission line ratios of small radio sources are in agreement with theoretical 
shock ionisation predictions and their velocity profiles are distorted. Together with 
. the other emission line properties this indicates that shocks associated with the ra- 

' dio source dominate the kinematics and ionisation of the emission line gas during 

I the period that the radio source is expanding through the interstellar medium. Gas 

^\ . clouds are accelerated by the shocks, giving rise to the irregular velocity structures 

observed, whilst shock compression of emission line gas clouds and the presence of the 
ionising photons associated with the shocks combine to lower the ionisation state of 
the emission line gas. By contrast, in larger sources the shock fronts have passed well 
O [ beyond the emission line regions; the emission line gas of these larger radio sources 

has much more settled kinematical properties, indicative of rotation, and emission 
^ ' line ratios consistent with the dominant source of ionising photons being the active 

' galactic nucleus. 

This strong evolution with radio size of the emission line gas properties of pow- 
i^j erful radio galaxies mirrors the radio size evolution seen in the nature of the optical- 

rS I ultraviolet continuum emission of these sources, implying that the continuum align- 

■ ment effect is likely to be related to the same radio source shocks. 

Key words: Galaxies: active — Galaxies: interstellar medium — Radio continuum: 
galaxies — Shock waves 
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1 INTRODUCTION 

Powerful high redshift {z > 1) radio galaxies display a num- 
ber of remarkable characteristics. Their near infrared emis- 
sion shows them to be amongst the most massive galaxies 
known in the early Universe and to have radial light profiles 
following de Vaucouleurs law, i ndicating that th ey are fully 
formed giant elliptical galaxies ( ^est et al. 1998 ). At optical 
and rest-frame ultraviolet (UV) wavelengths, however, the 
galaxies have very irregular morphologies, frequently show- 
ing a strong e xcess of emission aligned along the axis of th e 
radio source ( McCarthy et al. 1987] ; Chambers et al. 1987). 



Email: pbest@strw.leidenuniv.nl 



Their emission line properties are equally spectacular; lu- 
minous emission line regions surround the radio galaxies, 
extending for several tens of kiloparsecs or more (e.g. Mc- 
Carthy et al. 1995), with velocity shears up to a few hun- 
dred kms~^ and line widths as high as 1500 km s~^ (e.g. 
McCarthy et al. 1996). The emission line regions are charac- 
terised by high ionisation spectra including strong emission 
from species such as NeV and CIV. The origin of this lu- 
minous emission line gas, its kinematics and ionisation, and 
their connection to the radio source phenomenon remain im- 
portant astrophysical questions. 

Over the past few years we have been carrying out a 
detailed investigation of a sample of 28 radio g alaxies with 
redshifts z ~ 1 from the revised SCR catalogue (Laing et al 
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1983), using optical imaging with the Hubble Space Tele- 
scope (HST), high resolution radio interferometry with the 
Very Lar ge Array (VLA) and near infrared imaging with 



UKI RT ([Longair et al. 1995 



Best et al. 



1996| ; 



Best et al 



1997; Best et al. 1998). Here, and in an accompanying paper 



(Best et al. 1999; hereafter Paper 1), results are presented 
from a deep spectroscopic campaign using the William Her- 
schel Telescope (WHT) on 14 of these radio galaxies, to 
study in detail the emission line gas. The reader is referred 
to Paper 1 for details of the sample selection, data reduc- 
tion, the reduced one dimensional spectra and tabulated line 
fluxes, and the distributions of the intensity, velocity and 
line width of the emission line g function of position 

along the slit for each galaxy. 

The current paper is concerned with investigating the 
galaxy to galaxy variations in the ionisation and kinemat- 
ics of the emission line gas, and comparing these variations 
with the radio and optical properties. The layout is as fol- 
lows. In Section 2 the photoionisation and shock ionisation 
models are discussed and their predictions are compared 
with the observed emission line ratios of the galaxies on an 
emission line diagnostic diagram. The kinematics of the gas 
are studied in Section 3. In Section 4 the results are com- 
pared with those of a low redshift sample of galaxies and a 
scenario is proposed to explain the observed properties of 
the emission line gas at low and high redshifts. The spec- 
troscopic properties of the distant galaxies are compared 
to their optical properties, and the implications of these 
results for the continuum alignment effect are discussed. 
Conclusions are drawn in Section 5. Throughout the pa- 
per, values for the cosmological parameters of = 1 and 
Hq = 50kms"^Mpc-i are assumed. 



ising photons. The presence of such an obscured quasar nu- 
cleus is also indicated by the detection of spatially extended 
polarised emission and broad permitted lines in polarised 
light, due to scattering of the AGN light by electrons or 
dust (see Antonucci 1993 for a review). 

Photoionisation is not th e only story, however. As re- 
viewed by Binette et al. ( 1996 ) , simple photoionisation mod- 



els fail to reproduce some important features of the emission 
line spectra of the narrow line regions of active galaxies. 
In particular, the strengths of many high excitation lines 
(e.g. [NeV] 3426, CIV 1549, and high ionisation Fe lines) 
are under-predicted by factors as large as 10, the electronic 
temperatures derived by simple photoionisation models are 
too low when compared with those inferred from the line 
ratio [OIII] 4363 / [OIII] 5007, and photoionisation models 
alone cannot reproduce the large observed scatter in the 
Hell 4686 /H/3 ratio. Moreover, a significant fraction of ra- 
dio galaxies show indications of interactions between the 
radio jets and the surrounding emission line gas, with the 
radio source shocks determining the morphology and kine- 
matics of the gas. Detailed studies of individual sources (e.g. 
PKS 2250-41, z = 0.31, Clark et al. 1997, Villar-Marti'n 
et al. 1999; 3C171, z = 0.24, Clark et al. 1998) have shown 
that in some regions of the source the shocks can also dom- 
inate the ionisation; for example, minima in the ionisation 
state are observed coincident with the radio hotspots, and 
an anticorrelation is found between the ionisation state of 
the extended gas and its (jet shock broadened) line width. 

At high redshifts, z > 0.6, interactions between the ra- 
dio jets and the gas are readily apparent from the kinematics 
of the ionised gas (see Section ^). What has not been clear, 
however, is to what extent the shocks play a role in the 
ionisation of these high redshift sources. 



2 THE IONISATION OF THE EXTENDED 
EMISSION LINE REGIONS 



Robinson et al. (1987) showed that for most low redshift 
{z < 0.1) radio galaxies the emission line spectrum can 
be explained adequately if it is assumed that the gas is 
photoionised by a power-law emission source such as that 
provid ed by an active galactic nucleus (AGN). Baum et al. 
(1992) obtained similar results for a sample of radio galax- 
ies out to redshift 0.2, although they noted that photoion- 
isation models could not reproduce [Nil] 6584 /Ha ratios 
as high as were observed for some galaxies; they suggested 
that localised sources of heating and ionisation, for exam- 
ple shocks or the UV continuum of surrounding hot gas (e.g. 
Heckman e t al. 1 989), may play a role in some radio sources. 
McCarthy (1993) constructed a composite spectrum from a 



large sample of radio galaxies with redshifts 0.1 < z < 3; he 
showed that the emission line spectra of these more distant 
(more radio powerful) sources are also co nsisten t with being 
photoionised, but as Villar-Martin et al. (1997) argued, this 



composite spectrum is dominated by the few most highly 
ionised galaxies and so this result does not necessarily apply 
to the population as a whole. This photoionisation mecha- 
nism is in complete agreement with the currently popular 
orientation-based unification schemes of radio galaxies and 
radio loud quasars (e.g. Barthel 1989) in which all radio 
galaxies should host a powerful obscured active galactic nu- 
cleus, supplying a large flux of anisotropically emitted ion- 



2.1 The GUI] 2326 / GII] 1909 vs 

[Nelll] 3869 / [NeV] 3426 diagram 

Line r atio diagnostic diagrams, pioneered by Baldwin et al. 
(1981), provide a powerful tool for investigating the ionisa- 
tion mechanism of emission line gas. They have been widely 
used to distinguish the extended emission line regions of ac- 
tive nuclei from HII regions and planetary nebulae, and in 
recent years also between shock and photoionisation models 
for AGN. Standard emission line diagnostics at rest-frame 
optical wavelengths are shifted in to the near-infrared wave- 
bands for redshifts z > 1, and so cannot easily be used 
for high redshift radio galaxies. In the past couple of years, 
however, new diagnostic diagrams have been constructed for 
rest- frame UV emission lines redshifted into the o ptical win- 
dow dViUar-Martin et al. 19"97| ; [Allen et al. 1998| ). 

The emission line pairs CHI] 2326 & GII] 1909 and 
[Nelll] 3869 & [NeV] 3426 are well-suited for use as hue 
ratio diagnostics for distant radio galaxies for a number of 
reasons. All four of these fairly high excitation lines are 
relatively strong in AGN spectra, and so their fluxes can 
be determined with sufficient accuracy even for high red- 
shift radio galaxies. The two lines in each pair involve the 
same element, and therefore there is no dependence of the 
line ratios on metallicity or abundances; they are also close 
in wavelength, and so the effects of differential reddening 
or calibration errors are minimised. Perhaps most impor- 
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tantly, the predictions of shock and photoionisation models 
for these line ratios are significantly different. 

The line ratio Clll] 2326 / Cll] f 909 was determined for 
13 of the 14 sources in the sample from the data presented in 
Paper 1; the values of this ratio are given in Table |l] together 
with many other derived quantities of the radio galaxies. For 
the lowest redshift source, 3C340, CII] 1909 is not redshifted 
to a sufficiently high wavelength to be observed. The ratio 
[Nelll] 3869 / [NeV] 3426 is only available from the Paper 1 
data for the source 3C324, but for nine of the other sources 
in the sample values are available from the literature; these 
also are compiled in Table |l|. The two line ratios are plotted 
against each other on Figure |l|, where they are compared to 
the theoretical predictions of shock and photoionisation as 
discussed in the following subsections. 



2.2 Photoionisation models 

The theoretical line ratios of GUI] 2326 / CII] 1909 and 
[Nelll] 3869 /[NeV] 3426 for phot oionised gas were taken 
from the work of Allen et al. (1998) who calculated the pre- 
dicted line ratios for a number of emission li nes using the 
MAPPINGS II code ( [Sutherland et al. 1993| ). Allen et al. 
considered a planar slab of gas illuminated by a power-law 
spectrum of ionising radiation, and calculated the emission 
line ratios for a wide range of conditions: for two differ- 
ent spectral indices of the input spectrum [Fi, oc with 
Q = — 1 and a = —1.4), and two different densities of cloud 
(rie = 100 and 1000 cm~^), the ionisation parameter CTj^was 
allowed to vary in the range 10"* < [/ < 1. A high energy 
cutoff was applied to the ionising spectrum at 1.36 keV to 
avoid over-producing the intensity of the soft X-rays. The 
models are ionisation bounded and correspond to a range in 
cloud sizes from 0.003 to 32 parsec. The resultant line ratio 
sequences are shown on Figure |l| 

It is beyond the scope of this paper to provide a more 
detailed description of the photoionisation models, or of the 
other theoretical models considered next. For more complete 
discussions of these modelling techniques the reader is re- 
ferred to the papers from which the theoretical line ratio s 
have been drawn, in this case the work of Allen et al. (1998). 



2.3 Photoionisation including matter bounded 
clouds 

To avoid the shortcomings of simple photoionisation i nodel s 
discussed at the beginning of Section |2|, Binette et al. (199(: ) 
considered photoionisation of a composite population con- 
taining both optically thin (matter bounded; MB) and op- 
tically thick (ionisation bounded; IB) clouds (cf. Viegas and 
Prieto 1992). In their models, all of the photoionising radi- 
ation passes initially through the MB clouds which absorb 
a fraction (-Fmb ~ 40%) of the impinging ionising photons 
and produce the majority of the high ionisation lines in the 
spectrum. The radiation which isn't absorbed then strikes 



t The ionisation parameter U is defined as the ratio of the number 
density of ionising photons striking the cloud to the gas density 
(njj) at the front face of the cloud [U = (cnji)"^ / (F^du) /hu], 
where c is the speed of light and vo is the ionisation potential of 
hydrogen. 



the population of IB clouds; this radiation has already been 
filtered by the MB clouds as a result of which the IB clouds 
give rise to predominantly low and intermediate excitation 
lines. According to these models, the variation in the emis- 
sion line ratios from galaxy to galaxy has its origin in the 
variation of the ratio (hereafter Am/i) of the solid angle from 
the photoionising source subtended by MB clouds relative 
to that of IB clouds. A larger value of Am/i corresponds to a 
larger weight given to the MB clouds and hence a higher ex- 
citation spectrum. Since the model states that all of the ion- 
ising radiation striking the IB clouds must first have passed 
through the MB clouds, the ratio Ayi/i strictly cannot be 
below unity. 

Binette et al. consider two physical situations which 
might produce such a composite cloud population (some 
combination of the two would also be possible): (i) the MB 
'clouds' may be optically thin shells surrounding a denser 
IB core of a cloud; (ii) the MB clouds could be a separate 
population of clouds which lie close to the ionising source. 
In the second case, the MB clouds would have to have a 
covering factor of unity in order that all lines of sight to the 
IB clouds pass through a MB cloud; this is possible for very 
small clouds, consistent with the fact that they would be op- 
tically thin. Note that in the case of the MB clouds forming 
a separate cloud population, if some fraction of them are ob- 
scured from the observer, for example by the same material 
that obscures the active nucleus itself, this may give rise to 
an apparent Am/i < 1- 

Binette et al. tabulated the line ratios in the two cloud 
populations for a single set of parameters, chosen to give a 
good match to Seyfert spectra. They adopted a power-law 
spectrum with a spectral index of a = —1.3 {F,, oc u"), ion- 
isation parameters of [/mb = 0.04 and Um = 5.2 x 10~^, a 
density in the MB clouds of 50 cm^'^, and absorbed ionising 
photon fractions of -Fmb = 0.4 and Fib = 0.97 (see their pa- 
per for a more detailed discussion of these quantities). From 
these data, a sequence of line ratios of GUI] 2326 / GII] 1909 
and [Nelll] 3869 / [NeV] 3426 have been calculated allowing 
the quantity Am/i to vary in the range 0.001 < j4m/i < 100; 
this sequence is shown on Figure |l|. 

2.4 Shock ionisation models 

Fast radiative shocks are a powerful source of ionising pho- 
tons which can have a profound infiuence upon the tem- 
perature and ionisation properties of the gas in the post- 
shock region. An overview of shoc k ion isation models is pro- 
vided by Dopita & Sutherland (1996). The two most im- 



portant parameters for controlling the post-shock emission 
line spectrum are the velocity of the shock and the ratio 
B I ypn, where B is the pre-shock transverse magnetic field 
and n is the pre-shock number density of the emission line 
clouds. The latter ratio controls the density, and hence the 
effective ionisation parameter, of the post-shock gas since 
at high shock velocities the transverse magnetic field limits 
the compression caused by the shock through a balance be- 
tween the magnetic pressure of the cloud (oc B^) and the 
ram pressure of the shock (oc n) 



Dopita and Sutherland ( 1996 ) calculated the emission 
line ratios expected from gas ionised by the photons pro- 
duced in shocks for a range of physical conditions: the veloc- 
ity of the shock through the emission line clouds was allowed 
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Table 1. lonisation and kinematic properties of the radio galaxies, as calculated from the data presented in Paper 1. Column 1 gives the 
radio source name and column 2 its redshift. The projected linear size of the radio source (from Best et al. 1997) is given in column 3. The 
emission line ratios of CIII] 2326 / CII] 1909 and [Nelll] 3869 / [NeV] 3426, together with their uncertainties, are given in columns 4 to 7. 
The integrated [Oil] 3727 flux density is given in column 8 and its equivalent width in the rest-frame of the source is in column 9. The 
projected linear extent of the emission line region along the direction of the slit is in column 10. Column 11 gives the maximum FWHM 
of the [Oil] 3727 emission line in the spectrum, column 12 the range in relative velocities of this line seen along the slit, and column 13 
the number of distinct velocity components, as described in the text. Column 14 gives the extent of the aligned optical emission from 
the HST images of Best et al. (1997). 
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16 
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0.9 


1.2 
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0.3 
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78 
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800 
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24 
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1.079 


624 


8.4 


1.5 


0.9 


0.2 


0.76 


111 


32 


725 


50 
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12 


3C368 


1.132 


73 


1.1 


0.2 


3.9 


0.5 


5.87 
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3±1 
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0.47 
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900 
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Notes: 

[a] Values of the [Nelll] 3869 / [NeV] 3426 ratio are taken from: 3C22 — Rawlings et al. (1995); 3C217 & 3C340 — Spinrad, private 
communication; 3C265 & 3C352 — Smith et al. (1979); 3C280 — Spinrad (1982); 3C324 — Paper 1; 3C356 — Lacy & RawUngs 
(1994); 3C368 — Stockton et al. (1996); 3C441 — Lacy et al. (1998). 

[b] For 3C324 this corresponds to the two kinematically distinct components at +400 and — 400kms~^. 

[c] The data for 3C356 are for more northerly galaxy. 

[d] For 3C368 the determined equivalent width is a lower limit due to the contribution of the M-star (Hammer et al 1991) to the 
continuum level. 



to vary from 150 to 500kms~^, and the 'magnetic param- 
eter' was varied in the range < B/^Jn < 4/iGcm~^'^, 
which spans the expected range of values (see their paper 
for more details). These authors also emphasised the impor- 
tance of photons produced by the shock diffusing upstream 
and ionising the pre-shock gas. This may give rise to exten- 
sive precursor emission line regions, with different spectral 
characteristics to the compressed shocked gas. They there- 
fore calculated the emission line spectra predicted for these 
precursor regions for a range of shock velocities from 200 
to 500kms~^. In distant radio galaxies such as the ones 
studied in this paper, the spatial resolution is insufflcient 
to distinguish between the precursor and post-shock emis- 
sion regions, and a combined spectrum of the two would be 
observed. 

Using the data tables provided by Dopita & Suther- 
land, the emission line ratios of CIII] 2326 / CII] 1909 and 
[Nelll] 3869 / [NeV] 3426 were calculated both for simple 
shock models and for shock models including a precursor 
region. These theoretical ratios are shown on Figure |l|. 



are available for both lines clump into two groups of sources; 
3C217, 3C324, 3C352 and 3C368 are grouped together in the 
region corresponding to the shock ionisation models, while 
3C22, 3C265, 3C280, 3C340 and 3C356 are all close to the 
photoionisation predictions. Interestingly, the four sources in 
the first region all have projected radio linear sizes smaller 
than 115 kpc and the five sources in the second group all 
have sizes larger than this value. This result is more appar- 
ent in Figure | where the ratio of CIII] 2326 / CII] 19080 is 
plotted against the linear size of the radio source; the two 
parameters are correlated at the 98.5% significance level in 
a Spearmann-Rank test. 

For all nine sources, the photoionisation models of Bi- 
nette et al. (199C) including matter-bounded clouds provide 
an acceptable fit to the data. There is, however, a prob- 
lem with this model. As discussed in Section 2.3, a value of 



Am/1 < 1 is only possible if the MB clouds form a separate 
population of clouds, some of which are obscured from the 
observer. The four small radio sources would have to have a 
value Am /I ~ 0.1, implying that over 90% of their MB clouds 
are obscured, while there is no requirement for obscuration 



2.5 Interpreting the line diagnostic diagram 

Figure ^ clearly demonstrates that the ionisation states of 
radio galaxies, even within a tightly defined sample such 
as the one studied here, show considerable variations: the 
CIII] 2326 / CII] 1909 ratio differs by nearly a factor of ten 
between 3C368 and 3C340. The nine sources for which ratios 



t This ratio is preferred to the [Nelll] 3869 / [NeV] 3426 ratio 
here, and in later figures, because data are available for more of 
the galaxies and because the ratios are all drawn from the homo- 
geneous data set presented in Paper 1. The Neon ratio provides 
similar results, as is apparent from the strong inverse correlation 
between the two ratios seen on Figure hi 
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Figure 1. An omission line diagnostic plot for the SCR radio galaxies, compared with theoretical 
correspond to simple photoionisation models (a = —1.0 and a = —1.4), as described in Se ctio n 



sequence for photoionisation models including matt er bo unded clouds described in Section 2.3 



predictions. The upper shaded regions 
2l The dashed line corresponds to the 



The lower shaded region covers the ratios 
predicted by the shock models desc ribed in Section 2.4; the unshaded region just above this corresponds to the shock models including 
a precursor region (see Section 2. A). The five galaxies plotted towards the ed ge o f the diagram have no data available for one of their 
emission line ratios. An interpretation of this diagram can be found in Section 2.5. 



of the larger radio sources (Am/i ~ 1). The amount of ob- 
scuration would therefore have to depend upon the size of 
the radio source. 

Three factors determine the projected linear size of a 
radio source: the orientation of the source with respect to 
the line of sight, the age of the radio source and the advance 
rate of the hotspots. The first of these cannot be responsible 
since sources orientated more towards the line of sight, hence 
appearing smaller, will have less, not more, obscuration to- 
wards their central regions (cf. orientation-based schemes of 
radio galaxies and radio loud quasars) . Although the second 
option cannot be excluded, it seems improbable that the ob- 
scuration of MB clouds in the central regions will decrease by 
an order of magnitude during the short timescale of a radio 
source lifetime (a few x 10^ years) without the same process 
either destroying the clouds themselves or having other con- 
sequences, for example for the visibility of the broad-line re- 
gions. The third possibility, that there may be a connection 
between a higher obscuration of the central regions of small 
radio sources and a slower advance rate of their hotspots, has 
parallels with the suggestion that compact symmetric radio 
sources are small because they are confined by a dense (ob- 
scuring) surrounding medium (e.g. van Breugel et al 1984). 
However, recent investigations of compact radio sources us- 



ing VLBI techniques have derived hotspot advance veloc- 
ities of significant fractions of the speed of hght (~ 0.2c, 
e.g. Owsianik & Conway 1998, Owsianik et al 1998), sup- 
porting a youth rather than confinement scenario for these 
sources. This indicates that there is no strong connection 
between hotspot advance speed and obscuration (density) 
of the central regions of radio sources. 

It seems unlikely, therefore, that the value of Am/i 
should be strongly dependent upon the radio source size. 
Although this possibility cannot categorically be excluded, 
our preferred interpretation of the line diagnostic diagram 
is that for four sources, for which the extent of the emission 
line region is within a factor of two of the radio source size, 
the ionisation is dominated by shocks, whilst for the other 
five sources photoionisation dominates. This interpretation 
will be supported by the discussion of the kinematics in the 
following section; the matter-bounded cloud photoionisation 
model would provide no clear explanation of the variations 
seen in the kinematical properties. 

The uncertainties in the emission line ratios for each 
individual galaxy are too large to pin down any parame- 
ters of the ionisation accurately, but one feature is readily 
apparent. The five sources in the photoionisation region of 
the diagram are significantly more consistent with a flatter 
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spectral index (a ~ —1.0) for the power-law ionising contin- 
uum than with the steeper one (a ~ —1.4) typ ically adopted 
for low redshift sources. Villar-Martin et al. (1997) found a 
similar result analysing the rest-frame UV emission lines of 
radio galaxies with redshifts z > 1.7, suggesting that this 
may be a general feature of high redshift AGN. 



3 MORPHOLOGICAL AND KINEMATICAL 
PROPERTIES OF THE EMISSION LINE GAS 

The emission line gas surrounding low redshift radio galaxies 
shows velocity shears within the galaxies of between 50 and 
500kms~^, and (deconvolved) full width at half maxima 
(FWHM) of t he emission lines typica l ly in the range 20 



to 600kms~^( |Tadhunter et al. 1989bt |Baum et al. 1992[ ) 
In many cases the kinematics are consistent with a gravita- 
tional origin. At high redshifts the kinematics can be much 
more extreme, with velocity dispersions often in excess of 
1000 kms"^ (McCarthy et al. 1996, Paper 1) and compo- 
nents offset by several hundreds of kms~^ with respect to 
bulk of the gas (Tadhunter 1991, Paper 1). These remark- 
able kinematics are inconsistent with gravitational origins 
(cf. Tadhunter 1991). In this section the variation in the 
kinematics is compared with other properties of the radio 
source to investigate their origin. 

In Table ^ a number of the parameters of the emis- 
sion line properties of the gas are provided. The integrated 
[Oil] 3727 emission line intensity and the rest-frame equiv- 
alent width of this emission line are as calculated in Paper 
1. The projected linear size of the emission line region along 
the slit direction was determined from the extent of the lo- 
cations at which fits to the [Oil] 3727 emission line profile 
were obtained in Figures 2 to 15(d) of Paper 1 (excluding 
the detached emission line systems for 3C356 and 3C441). 
The range in relative velocities was calculated, to the nearest 
25kms~^, from Figures 2 to 15(e) of Paper 1, considering 
the velocity separation between the most positive and most 
negative velocity components of the [Oil] 3727 emission line, 
excluding any data points with uncertainties greater than 
100 km s"^ The maximum value of the FWHM of the emis- 
sion line gas was determined from Figures 2 to 15(f) in Paper 
1, again excluding any locations with uncertainties greater 
than lOOkms "\ 

One further parameter was calculated, hereafter re- 
ferred to as the 'number of velocity components', A'v, to 
provide an indication of the smoothness of the velocity pro- 
file. A^v was defined as the number of single velocity gradient 
components (ie. straight lines) necessary to fit, within the 
errors, the velocity profiles along the slit direction (Figures 2 
to 15e of Paper 1; cf van Ojik et al. 1997). A galaxy whose 
mean motion is consistent with simple rotation will provide 
a single component fit; higher values of A^'v correspond to 
irregular motions. This analysis, being by its very nature 
somewhat subjective, was carried out separately by two of 
the authors and by a third independent scientist. For 11 
of the 14 galaxies a unanimous value of A'^v was obtained. 
The remaining three galaxies show more complicated pro- 
files and their classification is ambiguous: for 3C280, values 
of 3, 4 and 5 were obtained, and so a value of 4± 1 is adopted; 
for 3C324 the profile is very different from the other galax- 
ies, being composed of two kinematically distinct systems as 



discussed in Paper 1 — a value of 2 is used; values of 2, 3 
and 4 were assigned to 3C368, and so 3 ± 1 is adopted. The 
precise values of A^v for these galaxies are of less importance 
than the fact they they are clearly inconsistent with a value 
of 1. The values of A'v are compiled in Table |l]. 

A number of features are immediately apparent from 
Table 0. It is noteworthy that of the seven galaxies with 
projected radio sizes smaller than 150 kpc, six have values 
A'v > 1 with only one having A'v = 1, while six of the seven 
sources larger than this size have A'^v ~ 1 (see Figure |^ . A 
chi-squared test shows that the probability of this occurring 
by chance is below 1%. Small radio sources predominantly 
have emission line gas with distorted velocity profiles and 
the emission line gas of large radio sources has a velocity 
profile generally consistent with rotation. 

A similar result is found with the variation of the max- 
imum FWHM with radio size, shown in Figure ^. These two 
parameters are anti-correlated at greater than the 99% sig- 
nificance level (Spearmann Rank test), with the four sources 
lying in the 'shock' region of the line diagnostic diagram 
(Figure ^) having clearly the highest values. This latter point 
is made more clearly in Figure |^ where the FWHM of the 
[Oil] 3727 emission can be seen to be inversely correlated 
with the CHI] 2326 / CII] 1909 emission line ratio, at the 
98.5% significance level using a Spearmann Rank correla- 
tion test. The kinematical and ionisation properties of these 
galaxies are fundamentally connected. 

It is not only the kinematics of the gas that evolve with 
the radio source size, but also the physical extent and the lu- 
minosity of the line emission. Figure y shows the variation of 
the equivalent width of the [Oil] 3727 emission line with in- 
creasing size of the radio source. Although this correlation is 
less strong (96% significance in a Spearmann Rank test) , it is 
apparent that the small sources in the 'shock-dominated' re- 
gion of Figure ^ show enhanced [Oil] 3727 equivalent widths. 
A more accurate description of Figure ^ is not that there 
is an inverse correlation between the equivalent width of 
the [on] emission and radio size, but rather that at large 
(> 150 kpc) radio sizes the distribution of equivalent widths 
is fairly fiat, and at small sizes there is often a factor of 2 to 
3 excess emission relative to this level. 

This enhancement of the line equivalent widths of small 
radio sources with respect to large sources implies an even 
greater boosting of their line luminosities, for two reasons. 
First, the optical continuum emission of small radio sources 
is more lumino us th an that of large sources, as indicated 
by Best et al. (1996), decreasing the apparent increase in 



the emission line equivalent width. Second, the equivalent 
width is determined from the extracted 1-dimensional spec- 
trum from a spatial region along the slit of about 35 kpc 
(see Paper 1); the physical extent of the emission line re- 
gions of small radio sources is greater than that of large ra- 
dio sources, as shown in Figure ^. Excluding 3C265, which 
is an exceptional source in many ways (e.g. see discussion in 
Best et al. 1997), the emission line regions of radio sources 
with sizes > 200 kpc have total extents of up to about 50 kpc 
(25 kpc radius, if symmetrical). Smaller radio sources, how- 
ever, have emission line regions ranging from this size up to 
about 100 kpc, a size comparable to the extent of the radio 
source. In other words, line emission at distances from the 
AGN of 30 to 50 kpc generally is only seen at the stage of ra- 
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Figure 2. The correlation between the CIII] 2326 / CII] 1909 
emission line ratio and the projected hnear size of the radio 
source. The four sources lying in the 'shocks' region of the line 
diagnostic diagram, Figure are plotted as filled triangles and 
the remainder of the galaxies as asterisks. 
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Figure 5. A plot showing the direct connection between the 
ionisation state of the emission line g indicated by the 

CIII] 2326 / CII] 1909 line ratio, and its kinematics in terms of 
the emission line maximum FWHM. Symbols as in Figure |2| 
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Figure 3. A histogram of the radio size distribution of the 
sources, separated into sources with smooth velocity profiles 
(A^v = 1, unshaded) and those whose profiles are irregular 
(Nv > 1, shaded). 
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Figure 6. The decrease in the equivalent width of the [Oil] 3727 
emission line with increasing size of the radio source. Symbols as 
in Figure jij 3C368 is plotted as a lower limit due to the contri- 
bution of the M-star to its continuum level. 
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Figure 4. The inverse correlation between the maximum FWHM 
of the [on] 3727 emission line and the projected linear size of the 
radio source. Symbols as in Figure eI 
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Figure 7. The variation of the linear extent of the [Oil] 3727 
emission line region with the size of the radio source. Symbols as 
in Figure ^ 
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dio source evolution when the hotspots are passing, or have 
just passed, through this region. 



4 DISCUSSION 

A number of resuhs have been derived in the previous sec- 
tions and these are summarised here for clarity. 

• Radio sources with small linear sizes (< 120 kpc) have 
lower ionisation states than larger radio sources. Their emis- 
sion line ratios are in agreement with the theoretical predic- 
tions of shock ionisation models, whilst those of large radio 
sources are consistent with photoionisation. 

• There is a strong inverse correlation between the 
FWHM of the [Oil] 3727 emission and the size of the radio 
source. The four sources with 'shock-dominated' ionisation 
states have the highest FWHM. 

• Large radio sources often have smooth velocity profiles 
consistent with rotation, whilst those of small sources are 
more distorted. 

• The [Oil] 3727 emission line strength correlates in- 
versely with the radio source size. The four 'shock- 
dominated' sources have the highest integrated [Oil] 3727 
equivalent widths. 

• The physical extent of the line emitting regions is larger 
in smaller radio sources. 

Before discussing the interpretation of these correla- 
tions, first a comparison is made to see if such results also 
hold for low redshift radio galaxies. 



4.1 Comparison with low redshift radio galaxies 



Table 2. Selected properties of the low redshift sample of FRII 
radio sources studied by Baum et al. (1992). 



Baum et al. (1992) studied the ionisation and kinematics of 
a sample of 40 radio galaxies with redshifts z < 0.2. Their 
sample contained a large mixture of radio source types, in- 
cluding both Fanaroff & Riley (1974; hereafter FR) class I 
and II object^ as well as sources with intermediate struc- 
tures. Many differences are now known to exist between the 
FRI and FRII sources besides the large differences at ra- 
dio wavelengths, such as the luminosity and environments 



of their host galaxies (Hill & Lilly 1991; Baum et al. 1995 
Ledlow fc Owen 1996), the luminosity of their emission line 
gas (Zirbel fc Baum 1995), differences in the dust proper- 



ties ( do Koff et al. 1999| ), and possib ly a different mode o f 
accretion on to the central black hole (Reynolds et al. 1996). 

Baum et al. also found a significant difference in the 
host galaxy kinematics between the two radio source types. 
They classified the kinematics of the radio galaxies into three 
classes, 'rotators', 'calm non-rotators' and 'violent non- 
rotators'. They found that almost all of the FRII sources fell 
into the rotator or violent non-rotator classes; most of the 
FRI and intermediate type sources were calm non-rotators. 
All of the FRH's had strong emission lines with a relatively 
high ionisation parameter, whilst the FR I and intermediate 
class sources had much weaker emission lines of lower ionisa- 
tion, with the surrounding hot interstellar and intracluster 



3 FR I radio sources are edge-darkened sources of generally lower 
radio luminosity that the FRII sources; FRII's are characterised 
by bright hotspots towards the extremities of each lobe. 
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15 
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Notes: 

[a] R — rotator; CNR — calm non-rotator; VNR — violent 
non-rotator. 

[b] No error quoted as the mean value is taken from data at only 
one position. 

medium likely to play an important role in the ionisation, 
both through heat conduction and through ionisation by 
its ultraviolet and soft X-ray emission (see also Baum et al. 
1995, Zirbel et al. 1995). Given the large differences between 
FRI and FRII sources, to allow a direct comparison with 
the high redshift sample, attention here is restricted to only 
the FRII sources in their sample. Four FRH's in the sam- 
ple lie more southerly than declination —30° and for two 
of these accurate determinations of the radio size could not 
be found in the literature; to avoid introducing any biases 
by selecting only the well-studied sources, all four of these 
sources have been excluded from further consideration. 

The remaining sample of low redshift FR II radio galax- 
ies is listed in Table |^, along with the linear size of the radio 
source taken from the literat ure a nd the kinematic classifi- 
cation given by Baum et al. (1992). In Figure]^ a histogram 
of the linear sizes is presented, separating the rotator and 
non-rotator classes. It is clear that the non-rotator classes 
are associated with small radio sources, and the rotator class 
with larger sources, exactly as is found for the high redshift 
sample. The only exception to this rule is 3C305, which is 
a rotator with a small radio size: indeed, this is the small- 
est radio source in the sample (15 kpc), and it could be ar- 
gued that any shocks associated with the radio source have 
not yet passed through a significant proportion of the host 
galaxy, accounting for the lack of clear non-rotational kine- 
matics. Even including 3C305, the probability of the radio 
sizes of the rotator and non-rotator classes being drawn from 
the same parent samples is less than 0.5% (using a Mann- 
Whitney U-test). 

All of the FRH's in the sample have relatively high 
ionisation states, bu t diff erences are seen from galaxy to 
galaxy. Baum et al. (1990) present the line strengths of the 
[01] 6300.3, [Nil] 6548.1,6583.4, Ha and [SII] 6716.4,6730.8 
emission lines as a function of position for half of the sample 
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Figure 8. A histogram of radio sizes for the different kinematic 
classifications of the FR II radio galaxies in the low redshift Baum 
et al. sample. 
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Figure 9. A plot showing the anticorrelation between the radio 
size and the ionisation state for the FRIFs in the low redshift 
Baum et al. sample. The rotators are represented by the crossed 
circles and the non-rotator by the filled triangles. 



considered in their 1992 paper. Although these are all rel- 
atively low ionisation lines and therefore not the most sen- 
sitive to differences between shock and photoionisation, the 
[OI] / Hq ratio should be somewhat higher for shock ionised 
gas than for photoionised gas. An 'average' value of this 
emission line ratio has been calculated for each galaxy as 
the mean of t he rat ios at the various positions tabulated by 
Baum et al. ( |1990| ); these are given in Table the errors 
quoted representing the scatter in the ratio with location in 
the galaxy. In Figure ^ these ratios are plotted against radio 
size: a Spearmann rank test shows that this emission line 
ratio is anticorrelated with radio size at the 96% confidence 
level. 

The low redshift sample therefore provides similar re- 
sults to the high redshift sample. Large FRII radio sources 
have kinematics consistent with rotation and higher ioni- 
sation states than small radio sources, whose ionisation and 
kinematics show more evidence for the role of shocks. It is of 
note that the low and intermediate redshift sources for which 
individual studies have shown unambiguously that the kine- 
matics and ionisation are dominated by shocks are almost 
invariably cases in which the radio source is of comparable 



size to the extended emission line regions (e.g. Clark et al. 
1997,1998), naturally agreeing with this picture. 

One significant difference that remains between the low 
and high redshift samples is that the high redshift sources 
are more extreme in their emission line properties (lumi- 
nosities, line widths, etc) than those at low redshifts. The 
most important factor influencing this is the sharp increase 
of the radio power with redshift in the fiux-limited samples, 
with corresponding increases in both the fiux of ionising pho- 
tons from the AGN and the energy of the jet shocks. How- 
ever, Tadhunter et al. ( 1998 ) investigated the correlations 
of different emission line strengths with redshift and showed 
that this cannot be the only reason: the ionisation-sensitive 
[Oil] 3727 / [OIII] 5007 ratio does not decrease strongly with 
redshift as it should if the only difference between the low 
power, low redshift and the high power, high redshift objects 
was that the latter contained a more luminous photoionis- 
ing source. They concluded that a secondary effect such as 
an increase in the density of the intergalactic medium or 
an increase in the importance of jet-cloud interactions with 
redshift is also required. 

4.2 The role of shocks in small sources 

The results presented in the previous sections, coupled with 
the evidence that a similar situation is seen at low redshifts, 
lead naturally to a single scenario to explain all of the emis- 
sion line properties. 

For small radio sources the morphology, kinematics and 
ionisation properties of the emission line gas are dominated 
by the effects of the bow shock associated with the expan- 
sion of the radio source. As this bow shock passes through 
the interstellar and intergalactic medium (ISM & IGM), the 
inter-cloud gas is quickly accelerated to the velocity of the 
bow shock, but the warm emission line clouds are essen- 
tially bypassed by the shock front (e.g. Rees 1989, Begel- 
man & Cioffi 1989). The clouds are accelerated during the 
short time it takes the shock to pass the cloud by, due to the 
imbalance in the pressures between the pre-shock and post- 
shock gas on the front and back of the cloud. The velocity to 
which the clouds are accelerated in this way is easily shown 
to be independent of cloud size and to be well below 100 kms 
(e.g. Rees 1989). 

Much larger velocities are induced, however, if the ef- 
fect of ram-pressure acceleration by the shocked IGM gas 
(often referred to as entrainment) is considered. Behind the 
initial bow shock, the clouds find themselves in a shocked 
layer of IGM, moving outwards at speeds approaching that 
of the bow shock. The clouds will be accelerated within this 
medium until they pass across the contact discontinuity into 
the radio cocoon, where the pressure is the same as in the 
shocked layer of gas but the density is much lower, and they 
are no longer accelerated; there is essentially no mixing of 
the hot inter-cloud gas across this contact discontinuity (e.g. 
Norman et al. 1982). 

During the time AT for which a cloud is between the 
bow shock and the contact discontinuity, the momentum 
imparted to the cloud by the shocked IGM can be approxi- 
mated to first order as r^v^UgnipAT, where is the cloud 
size, Vs is the bow shock velocity, rig is the post-shock num- 
ber density of the inter-cloud gas and rup is the proton 
mass. The mass of the cloud is of order r^Ucmp, where ric is 
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the cloud number density, and the timescale AT is of order 
D/vs, where D is the distance between the bow-shock and 
the contact discontinuity; therefore, the velocity to which 
the cloud is accelerated is 
D tIe 



(1) 



In the radio source evolution models of Kaiser 



Alexander (1997), radio sources grow self-similarly and D 



is found to be about 3% of the distance between the AGN 



and the bow shock (Kaiser & Alexander 1999); for the radio 



source passing through the emission line region at radius 
~ 15kpc then, D ~ 0.5 kpc. Assuming a density ratio of 
Tie /uc ~ 10"* for pressure equilibrium between the clouds 
(T ~ 10" K) and the surrounding IGM (T ~ 10* K), and a 
shock velocity of Vg ~ 0.05c from typical hot-spot advance 
velocities (e.g. Liu, Pooley & Riley 1992), then a cloud of 
size Tc ~ 1 pc will be accelerated to 750 kms comparable 
to the velocities observed in the small radio sources (the ac- 
tual velocities may need to be slightly higher since the radio 
galaxies are believed to lie close to the plane of the sky) . The 
spread in projected cloud velocities from clouds of different 
sizes and from clouds accelerated through different regions 
of the bow shock will lead to the broad velocity dispersions. 
The acceleration of the emission line gas clouds by the radio 
bow shocks therefore explains the distorted velocity profiles 
and large line widths observed in small radio sources. 

It is interesting to note that the acquired cloud veloc- 
ities are proportional to the bow-shock velocity Vs- If the 
bow-shock velocity increases with radio power (redshift) , as 
has been suggested from spectral ageing measurements of 
hotspot advance velocities (e.g. Liu, Pooley & Riley 1992), 
this would explain why greater velocity widths are seen in 
high redshift sources than low redshift sources. 

The ionisation state of the large radio sources indicates 
that the dominant source of ionising photons in these sources 
is the AGN. Since the properties of the AGN are not ex- 
pected to change dramatically between small and large radio 
sources, the gas surrounding the small sources should receive 
a similar flux of photoionising radiation. The lower ionisa- 
tion state seen in the spectra of these galaxies arises in part 
due to compression of emission line gas clouds by the ra- 
dio source shocks, decreasing the ionisation parameter. The 
presence of extra (softer) ionising photons associated with 
the shocks further influences the ionisation state. Bicknell 
et al. (1997; see also Dopita and Sutherland 1996) have in- 
vestigated the emission line luminosity that can be generated 
by radio source shocks expanding through a single-phase 
ISM with a power-law density gradient, p(r) = po^r /ro)~^ ■ 
For 5 — 2, they show that the work done on the ISM by the 
expanding radio cocoon (PdV^) is approximately half of the 
energy supplied by the radio jet; if the shock is fully radia- 
tive then a signiflcant proportion of this energy is fed into 
emission line luminosity. The luminosity of the [OIII] 5007 
emission line can be estimated as 



Vio-iv V] 



PlA 



X 10*"^ ergs s ^, 



^([0™])--8-5VlO-y VlO-WHz 
where Pi. 4 is the monochromatic power of the radio source 
at 1.4 GHz, and kia is the conversion factor from the energy 
flux of the jet to the monochromatic radio power at 1.4 GHz, 
which BickneU et al. ( |l997[ ) estimate to be of order 10"^" ^ 
Adopting this value, taking the flux density of a typical z ~ 



1 SCR source at an observed frequency of 1.4 GHz to be 

2 Jy, and assuming that t he [Oil] 3727 / [Q UI] 5007 emission 
line flux ratio is ~ 0.5 ( [McCarthy 1993D , then for 5 = 2 
the observed [OH] 3727 emission line flux produced by the 
shocks is calculated to be /([OH]) ~ 3x 10~^^ ergss"^ cm"'^. 
Of this, probably between a third and a half (that is, 1 to 
1.5 X 10~^^ ergss"^ cm~^) will fall within the projected sky 
area from which the spectrum was extracted. This predicted 
emission line flux can be compared to the [Oil] 3727 emission 
line fluxes observed in the data (Table |l|), which lie in the 
range ~ 0.5 to 5 x 10~^^ ergs s~^ cm~^ with the smaller radio 
sources typically having the higher values (see also Figure ^ . 
These results are completely consistent with a small (factors 
of 2 to 5) boosting of the emission line luminosities of small 
sources due to the extra energy input from the shocks. 

Once the radio source shocks have passed beyond the 
emission line clouds, the shock-induced emission line lumi- 
nosity will fall. Under the simplest assumptions, once the 
jets pass beyond the confining ISM the pressure inside the 
cocoon will drain away, and the cocoon wall shocks will no 
longer be pressure driven (e.g. Dopita 1999). These shocks 
will pass into a momentum conserving phase; their velocity 



will decrease roughly as Vs 



and so since the shock lu- 



minosity per unit area scales as , the shock-induced emis- 
sion line luminosity will fall as r"'*. Although these assump- 
tions are oversimplified, taking no account of confinement by 
an intracluster medium for example, it is clear that once the 
shock fronts have now passed well beyond the emission line 
regions, the contribution of ionising photons produced by 
the shocks will decrease rapidly; this is in complete accord 
with the larger sources having photoionisation dominated 
emission line regions. 

4.3 The physical extent of the emission line gas 

The physical extent of the emission line region of each galaxy 
along the slit direction was provided in Table ^. For com- 
parison, the extent of the aligned optical (rest-frame ultra- 
violet) emission has also be en de termined from the HST 
observations of Best et al. (1997); using the HST image 



taken through the filter at a rest-frame wavelength of about 
4OOOA, the angular distance over which optical emission was 
observed at greater than three times the rms sky noise level 
of the image was measured for each galaxy, and the corre- 
sponding 'optical sizes' are given in Table 1. These values 
can fu rther be compared with these results of Best et al. 
( [1998| ), who showed from near-infrared imaging that, un- 



derlying the aligned emission, the radio sources are hosted 
by giant elliptical galaxies with characteristic radii of typi- 
cally 10 to 15 kpc. 

The extent of the optical aligned emission does not ex- 
ceed 25 kpc except in three cases: 3C247, 3C265 and 3C368. 
For 3C368, the HST 'continuum' image is actually domi- 
nated by a combination of line emission and the correspond- 
ingly lumi nous nebular continuum emission (see discussion 
in Section 4.5). The large extent of 3C247 is also likely to 



be predominantly line emission, since it arises from a diffuse 
halo of emission exactly tracking that seen i n a narrow-band 
[OH] 3727 image by McCarthy et al. ( |l99£| ). 

With the exception of 3C265 (which, as discussed in 
Paper 1, is an unusual source in many ways), it is therefore 
reasonable to say that the aligned continuum emission has 
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an extent of only a couple of characteristic radii, and so lies 
within the body of the host galaxy. The situation with the 
emission line gas is very different: this has a physical ex- 
tent which can exceed 100 kpc, with a mean extent of over 
60kpc. The emission line gas clearly extends well beyond 
the confines of the host galaxy. As was shown in Figure ^ 
there is also a difference in the physical extent of the line 
emitting regions between large and small radio sources, with 
line emission at radii of 30 to 50 kpc generally only seen in 
small radio sources. Unless there is an intrinsic difference be- 
tween the environments of the small and large radio sources, 
which seems unlikely given all of the correlations found, the 
emission line gas clouds must also be present out to radii 
> 30 kpc in large radio sources, but is not visible. 

Again, the role of shocks can be considered to explain 
this. At these radii, the flux of ionising photons from the 
active nucleus may be insufficient to produce an observable 
emission line luminosity. As the radio source shocks pass 
through these regions, however, the gas density will be in- 
creased and, as discussed above, a large source of local ionis- 
ing photons will become available, pushing up the emission 
line luminosity. Following the passage of the radio shocks 
and the consequent removal of the associated ionising pho- 
tons, this enhanced line emission will fade over timescales 
much shorter than the radio source lifetime. Thus, luminous 
line emission is only seen from the clouds at radii 30 to 
50 kpc at the time that the radio source shocks are passing 
through these regions. A direct consequence of this model 
is that for radio sources smaller than about 100 kpc a posi- 
tive correlation between radio source size and emission line 
region size should be observed, since line emission from the 
clouds at radii 30 to 50 kpc will not be seen until the radio 
source has advanced that far. Such a correlation has indeed 
been o bserved in the Ly-a emission of radio galaxies with 
z ^ 2 ( Ivan Ojik et al. 1997|) . 

An interesting test of the model presented here could be 
carried out by taking high spatial resolution long-slit spec- 
tra of a sample of radio galaxies with radio sizes smaller 
than the size of the emission line regions. The prediction 
is that within the region of the host galaxy occupied by 
the radio source, the radio source shocks will be important; 
the emission line ratios will be consistent with shock ioni- 
sation, and the gas kinematics will be distorted with broad 
velocity dispersions. Outside of this region, however, the gas 
clouds will not yet have been influenced by the radio source 
shocks and photoionisation should dominate. A study with 
a similar principle has been carried out on the radio source 
1243+036, a radio galaxy of radio size about 50 kpc at a 
redshift z = 3.6. Distorted Ly-a velocity structures with 
large velocity FWHM are seen within the radio source struc- 
ture, but Ly-Q emission also extends beyond that to at least 



75 k pc r adius in an apparent ly rot ating halo ( van Ojik et al _ 
199e|)7^illar-Marti'n et al. ( |l999[ ) have also found that the 
line emission of PKS 2250—41 {z — 0.308) is composed of 
distinct kinematic components; a low ionisation component 
with broad velocity width in the region of the radio source 
structure, and a narrower high ionisation component which 
extends beyond the radio lobe. Carrying out studies such as 
these for a large sample of radio sources is important because 
the velocity structures of the line emission in regions outside 
the radio shocks will directly show the initial motions of the 
emission line clouds and can be used to determine whether 



these clouds are simply material associated with the forma- 
tion of the galaxy which has been expelled into the IGM, or 
whether they have an external origin, brought in either by a 
galaxy merger or a cooling flow. It is difficult to distinguish 
between such scenarios in larger radio sources since infor- 
mation on the initial cloud velocities has been destroyed by 
the bow shock acceleration. 



4.4 Evolution of the velocity structures 

One significant issue remains to be explained in this picture, 
and that is how the velocity structures of the large radio 
sources are produced. The high gas velocities and velocity 
dispersions induced by the shocks in small radio sources are 
seen to evolve within the timescale of a radio source life- 
time, a few X 10^ years, such that the emission line clouds 
obtain an underlying velocity profile consistent with a rotat- 
ing halo, albeit with a still high velocity dispersion. Ques- 
tions that need to be considered are whether this truly is 
rotation that is being seen, over what timescale can the ex- 
treme shock-induced kinematics be damped down, and can 
a mean rotation profile be produced whilst the velocity dis- 
persion remains so high? 

Regarding the first question, given the single slit posi- 
tion and relatively low spatial resolution for the high red- 
shift radio galaxies, it cannot categorically be stated that 
the emission line profiles of large radio sources are rotation 
profiles. The data are consistent, for example, with outfiow 
along the radio axes, although in this case it is not clear why 
the velocity increases with radius (a structure more like that 
of 3C324 - see Paper 1 - might be expected) while the ve- 
locity dispersion decreases with increasing source size. At 
low redshifts, however, much higher spatial resolution stud- 
ies using multiple slit positions show clearly that the gas is 
in rotating structures (e.g. Baum et al. 1990). It therefore 
seems reasonable to assume that this may also be true at 
higher redshifts, and even if this is not the case, the questions 
noted above still need to be addressed for the low redshift 
radio sources. 

Three plausible mechanisms can be considered for the 
evolution in the velocities of the emission line clouds over 
the radio source lifetime. The first is that the emission line 
clouds settle back into stable orbits within the host galaxy 
through gravitational dynamics alone. The timescale for this 
process is of order a few crossing times of the clouds, where 
for clouds moving with velocity Vc ~ 500 km s ~^ at a radius 
r ~ 15 kpc in the galaxy, the crossing time is tc ~ 2r/v ~ 
6 X lO'^ years. This timescale is longer than the radio source 
lifetime, and so gravity alone cannot give rise to the observed 
evolution in the emission line structures. 

A second possibility concerns the deceleration of emis- 
sion line clouds moving with respect to the interstel- 
lar medium, due to ram-pressure arguments. This works 
through the same process as the acceleration argument dis- 
cussed in Section 



1.2 



As the emission line clouds move 
through the inter-cloud gas, those clouds moving with the 
largest velocities sweep up the greatest mass of inter-cloud 
gas, and so are decelerated most strongly. This process will 
decrease the width of the cloud velocity distribution. 

Simulations have been carried out, as detailed in Ap- 
pendix to investigate the timescale over which the mean 
velocity of an ensemble of emission line clouds (with an ini- 
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tial velocity distribution similar to that seen in small radio 
sources) evolves to that of the IGM in which the clouds 
are moving, and the timescale over which the dispersion of 
the velocity distribution is decreased. It is found that the 
peak of the velocity distribution evolves to that of the gas 
in which it is moving considerably more quickly than the 
velocity width decreases. Both timescales depends upon the 
typical cloud size and the ratio of the cloud density to that 
of the inter-cloud medium within the radio cocoon, and for 
reasonable assumptions the timescale for decrease of the ve- 
locity widths is comparable to the radio source lifetime (see 
Appendix ^ for details) . 

Therefore, if a population of emission line clouds were 
placed within the rotating ISM of a galaxy, a mean rotation 
profile for the emission line clouds could be recovered whilst 
the FWHM of the emission lines remained large, as is ob- 
served in the radio galaxies. The problem with this model, 
however, is that the radio bow shock sweeps up essentially 
all of the inter-cloud gas, with little mixing through the 
contact discontinuity (e.g. Norman et al. 1982). The radio 
cocoon is filled primarily with material supplied by the ra- 
dio jets, and so there is essentially no gas left following a 
rotation profile. Such gas would have to be resupplied to 
the ISM, for example by supernovae and stellar winds from 
stars in rotational orbits, but it is very unlikely that enough 
gas can be supplied in this manner. Alternatively, the co- 
coon material supplied by the radio jets would itself have to 
be in rotational motion, perhaps through angular momen- 
tum transfer from the rotating IGM to the radio source as 
the bow shocks advance. To summarise, although this mech- 
anism can decrease the velocity widths of the gas, it is not 
clear whether a rotation profile can be re-established quickly 
enough. 

The third possibility considers the evolution of the pop- 
ulation of radiating clouds through a combination of galaxy 
rotation and cloud shredding. Klein et al. (1994) showed 
that in the aftermath of a bow shock, emission line clouds 
may be susceptible to shredding due to growing Kelvin- 
Helmholtz and Rayleigh-Taylor instabilities on their sur- 
face. The clouds could be shredded over a timescale of a 



few 'cloud crushing times', tc 



X 



1/2, 



;/wb, where x is the 



density ratio of the cloud to the surrounding medium in- 
side the cocoon, is the post-shock cloud radius, and lib is 
t he ve locity of the bow shock through the IGM. Kaiser et al. 
(1999) considered such cloud disruption as a way to resupply 
material to the radio cocoon in order to explain how a sec- 
ondary hotspot can be formed in the newly discovered class 
of double-double radio galaxies (e.g. Schoenmakers et al. 
1999); they derived a value of tec ~ 5 x 10^(rc/pc) yrs. 

The cloud shredding time is shortest for the smallest 
emission line clouds; clouds smaller than about a parsec will 
be shredded on timescales shorter than the radio source life- 
time. These small clouds were the most rapidly accelerated 
(see Equation |l|) and so are responsible for producing much 
of the high velocity dispersions and distorted velocity struc- 
tures. If these high velocity small clouds are destroyed then 
the line emission will become dominated by the remaining 
more massive clouds, which were less accelerated by the ra- 
dio source shocks, have a lower velocity dispersion, and may 
still maintain the vestiges of a rotation profile. 

Equation ^ further shows that the velocity acquired by 
the warm clouds is proportional to the velocity of the bow 



shock. In directions perpendicular to the radio axis the bow 
shock velocity is lower by a factor of the aspect ratio of the 
cocoon (typically between about 1.5 and 6, e.g. Leahy et al. 
1989), and so the warm clouds in these directions will be less 
accelerated. In small radio sources these clouds will not be 
very luminous since they lie away from the strongest radio 
source shocks and outside of the cone of photoionising ra- 
diation from the partially obscured AGN; the emission will 
be dominated by the higher velocity clouds along the radio 
jet direction. Over a rotation timescale (~ lO'^ years), how- 
ever, these low velocity clouds may be brought within the 
ionisation cone of the AGN, become ionised, and contribute 
significantly to the emission line luminosity. Likewise, clouds 
in small radius orbits around the AGN will acquire lower ve- 
locities, since the distance between the bow shock and the 
contact discontinuity is less and so the period of acceleration 
is shorter. Thus the rotation profile may re-establish itself 
from the central regions of the galaxy outwards. 

By these two mechanisms of shredding and mixing of 
the cloud populations, the observed population of emission 
line clouds will evolve such that, in large radio sources, an 
increasing percentage of the emission will arise from clouds 
which were less accelerated by the bow shock and so the ro- 
tation profile will be gradually recovered. The cloud shred- 
ding model has a further advantage that if some fraction 
of the clouds are destroyed in large radio sources then the 
emission line luminosity will decrease with increasing radio 
size, as is observed. The one drawback of this model is that 
it is surprising that the distinction between radio sources 
showing rotation profiles and those with distorted profiles is 
so sharp. Note also that if this scenario is the correct one 
then the emission line clouds must lie in rotating orbits prior 
to the radio source activity, providing some information as 
to their origin. 

In conclusion, the observation of emission line clouds in 
rotating halos around large radio galaxies is not trivial to 
explain, given the large influence of the radio source bow 
shocks passing through the medium. Gravitational effects 
alone cannot be responsible for re-establishing rotation pro- 
files, but a combination of cloud shredding and cloud mix- 
ing, maybe with some help from ram-pressure deceleration, 
could reproduce the effect. 

4.5 Implications for the alignment effect 

In 1987, McCarthy et al. and Chambers et al. demonstrated 
that the optical-UV emission of radio galaxies with redshifts 
z > 0.6 has a strong tendency to be elongated and aligned 
along the direction of the radio sou rce. HST image s of a 
sample of 28 of these radio galaxies (Best et al. 1997) have 



demonstrated that the form of this so-called 'alignment ef- 
fect' varies strongly from galaxy to galaxy, and in particular 
appears to evolve with increasing size of the radio source 
(Best et al. 1996). Small radio sources show a number of 



intense blue knots tightly aligned along the direction of the 
radio jet, whilst larger sources generally have more diffuse 
optical-UV morphologies. Given the strong similarity be- 
tween this radio size evolution of the continuum alignment 
effect and the evolution of the emission line gas properties, it 
is instructive to examine the role of the radio source shocks 
and the emission line clouds in giving rise to continuum 
emission. 
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One direct connection is the nebular continuum emis- 



sion from the warm emission line gas clouds (Dickson et al 
1995), that is, free-free emission, free-bound recombination, 
two-photon continuum and the Balmer forest lines. The flux 
density of this emission is directly connected to the flux of 
the H/3 emission line. The very luminous line emission seen 
in the spectra of these powerful radio galaxies (e.g. Paper 



fraction of the excess optical-UV emission must be associ- 
ated with this mechanism. However, the lack of polarised 
emission from some sources (e.g. 3C368, van Breugel et al. 
1996; see also Tadhunter et al. 1997) dictates that this is not 
universal; even for 3C324 where the polarisation percentage 
is high, only a fraction < 30 — 50% of the optical- UV emis- 
sion is associated with the scattered component ( [Cimatti 



1) thus implies that nebular continuum emission is likely et al. 1996). A problem for scattering models is that, in the 



to make a significant contribution to their UV flux density. 
Indeed, 3C368 was one of the o riginal three radio galaxies 
studied by Dickson et al. ( 1995 ), and they found a nebular 
continuum contribution in the northern knots as high as 60% 
of the total continuum emission at rest-frame wavelengths 
just below the 3646A Balmer break (see also Stockton et al. 
1996). As can be seen from Table hi, 3C368 is a somewhat 
extreme case and the contribution for more typical galaxies 
will be somewhat lower, but still of great significance. In Sec- 
tion ^ it was shown that the luminosity of the emission lines 
correlated inversely with the size of the radio source (Fig- 
ure^; therefore, the strength of nebular continuum emission 
will decrease with increasing radio source size, and in small 
sources will be found predominantly along the radio jet trac- 
ing the strongest radio source shocks. This reflects exactly 
the observed evolution of the continuum alignment effect. 

A second alignment effect hypothesis involving the 
emission line clouds is that star formation is induced by 
the passage of the radio jet, due to the radio source shocks 
compressing gas clouds and pushing them over the Jean's 
limit (e.g. Rees 1989, Begelman and Cioffi 1989, De Young 
1989). It should be noted that it is the most massive clouds 
which would collapse to form the stars, and these are dis- 
tinct from the smallest clouds which are the most likely to be 
destroyed by the bow shock. In regions which might be star- 
forming, < 10® years behind the bow shock, the only clouds 
which will already have been destroyed by instabilities on 
their surface are those of size Vc < 0.1 pc (see Section 4.3); 
for a mean cloud density of 100 cm~^ this corresponds to a 
total cloud mass of less than lO~^Af0, not massive enough 
to have formed a star anyway. 



As discussed by Best et al. (1996), the jet-induced star 



formation mechanism can also account directly for the evolu- 
tion of the optical-UV morphology with radio size: the mass 
of stars required to pro duce the excess optica l -UV emission 



is o nly a few xlO Mq ( [Lilly fc Longair 1984| ; punlop et al 



198S |), y ell below 1% of the stellar mass of the galaxy, and 
since the starburst luminosity drops rapidly with age they 
become indistinguishable from the evolved star population 
over a timescale of a few x 10^ years. On the negative side, 
no direct evidence for young stars in these radio galaxies was 
found in our spectra (cf. 4C41.17 at higher redshift, z = 3.8; 
Dey et al. 1997), although the clearest features of young stel- 
lar populations fall outside the observed wavelength ranges. 

Another important continuum alignment model is scat- 
t ering of ligh t from a hidden quasar nucleus by electrons 
(Fabian 1989) or dust (e.g Tadhunter et al. 1989a, di Serego 
Alighieri et al. 1989). Strong support for this model comes 
from the observation that the optical emission of some dis- 
tant radio galaxies is polarised at the ~ 10% level with the 
electric vector oriented perpendicular to the radio axis (e.g. 
Cimatti et al. 1996 and references therein), a nd the detection 
of b road permitted lines i n polarised light ( Dey & Spinrad 
1996; C Imatti et al. 1991: ; Tran et al. 1996): clearly some 



simplest picture, a biconical emission region is expected for 
the scattered light, rather than the knotty strings of emis- 
sion observed to lie along the radio jet. However, in light of 
jet -shock models, this could be explained by extra scattering 
particles being made available along the radio jet axis, either 
as dust grains being produced in jet-induced star forming 
regions, or by radio source shocks disrupting optically thick 
clouds along the rad io jet direction and exposing previously 
hidden dust grains (Bremer et al. 1997). 

In conclusion, radio source shocks will play a key role 
in producing the observed morphology and radio size evolu- 
tion of the continuum alignment effect. Nebular contirmum 
emission will be enhanced in small radio sources, some gas 
clouds may be induced to collapse and form stars, and extra 
scattering particles associated either with any star formation 
or the disruption of gas clouds could enhance the scattered 
component. 



5 CONCLUSIONS 

The main conclusions of this work can be summarised as 
follows: 

• Small radio sources show a lower ionisation state than 
large radio sources. The emission line ratios of radio sources 
with linear sizes < 120 kpc are consistent with the gas being 
ionised by photons produced by the shocks associated with 
the radio source. The emission line luminosities of the small 
sources are boosted by a small factor (~ 2 — 5) relative to 
large sources, in accord with them receiving an extra source 
of ionising photons from the shock. 

• Small radio sources have very distorted velocity profiles, 
large velocity widths, and emission line regions covering a 
larger spatial extent than those of large sources; the lat- 
ter have much smoother velocity profiles which appear to 
be dominated by gravitation. These properties are fully ex- 
plained in terms of the passage of the shocks associated with 
the radio source. 

• A strong correlation is found between the ionisation 
state of the gas and its kinematical properties, indicating 
that the two are fundamentally connected. 

• These correlations, originally derived for the sample of 
redshift one radio galaxies studied in Paper 1, are shown also 
to hold for a sample of FRII radio galaxies with redshifts 
z < 0.2. 

• The similarity of the evolution of the emission line gas 
properties with radio size to that of the continuum alignment 
effect makes a strong case for the continuum alignment effect 
also having a large dependence upon radio source shocks. 

• The continuum alignment effect is generally confined 
to within the extent of the host galaxy, but line emission is 
observed over a considerably larger spatial extent. 
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APPENDIX A: DECELERATION OF EMISSION 
LINE CLOUDS MOVING THROUGH THE IGM 

Consider a spherical cloud of emission line gas with number 
density ric and radius rc travelling at velocity Vc through gas 
of number density rig and velocity iig. In a time dt a mass 
of gas of approximately ■Kr^ngTnp{vc — Vg)dt, where rrip is 
the proton mass, is displaced by the cloud and accelerated 
from velocity Vg to velocity Vc. The momentum of the cloud 
is correspondingly decreased: 

— nr^ricmpdvc = nr'^ngnip^Vc — Vg)^dt 

Defining tg as tg = 4ncro/3ng, where ro is a typical 
cloud radius, then 



dvc 



dt 



Tc/ro to 



Using this equation it is possible to follow the evolution 
of an ensemble of such emission line clouds. For simplic- 
ity the distribution of emission line cloud radii was chosen 
to be flat in logarithm space over a factor of 1000 range 
centred on ro, that is, P(log(rc)) is constant in the range 
— 1.5 < log(rc/ro) < 1.5, and outside that range. The 
initial velocity distribution of the clouds was set to follow 
a Gaussian distribution with a mean velocity of zero and 
a FWEIM of 1000 km s~^, chosen to represent the velocity 
dispersion observed in small radio sources. A Monte Carlo 
simulation was then used to follow the evolution of the ve- 
locity distribution of the cloud population in gas moving 
with velocity -1-200 km s~^, typical of the relative velocity 
offsets seen in the galaxy profiles; the results are shown in 
Figure 



Al 



It can be seen that the peak of the velocity distribu- 
tion of the cloud population evolves rapidly to that of the 
gas in which it is moving; the width of the velocity distribu- 
tion becomes progressively narrower but over a much longer 
timescale. The resulting velocity distribution is no longer 
Gaussian, but can be approximated as a Gaussian distribu- 
tion plus extended broad wings (the slight dip at 200 km s ~^ 
for the first plotted time interval should be ignored; it arises 
only due to the simplicity of the model in which a cloud 
whose initial velocity is close to 200kms~^ will decelerate 
very slowly). 

The timescale over which the FWHM of the emission 
line clouds decreases to that observed in large radio sources 
(a few hundred kms can be used to test the plausibility 
of this model. This time interval is T « 5 x 10~^to (the 
solid line on Figure Al), where to = 4ncro/3ng. Taking T ~ 
3 X 10^ yr as an appropriate age for radio sources a few 



T-5x10"'to 

T-2x10"% 

T=5,<10-=t„ 




Velocity (km/s) 

Figure Al. Modelling the evolution of the emission line cloud 
velocity distribution within gas moving at 200 km s ~ ^ . 



hundred kpc in size, and assuming pre-shock density ratio 
n-c/rig ~ IC* with the IGM density decreased a further factor 
A by the bow shock, this gives (ro/pc)A « 5. For A ~ 40 
(e.g. Clarke & Burns 1991) the median cloud size would be 
about about 0.15 parsec; although small, this is certainly 
plausible given the simplicity of the assumptions. 
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